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Influence of the heterointerface sharpness on exciton recombination dynamics in an 
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The dynamics of exciton recombination in an ensemble of indirect band-gap (In,Al)As/AlAs 
quantum dots with type-I band alignment is studied. The lifetime of confined excitons which are 
indirect in momentum-space is mainly influenced by the sharpness of the heterointerface between the 
(In,Al)As quantum dot and the AlAs barrier matrix. Time-resolved photoluminescence experiments 
and theoretical model calculations reveal a strong dependence of the exciton lifetime on the thickness 
of the interface diffusion layer. The lifetime of excitons with a particular optical transition energy 
varies because this energy is obtained for quantum dots differing in size, shape and composition. The 
different exciton lifetimes, which result in photoluminescence with non-exponential decay obeying 
a power-law function, can be described by a phenomenological distribution function G{t), which 
allows one to explain the photoluminescence decay with one fitting parameter only. 

PACS numbers: 78.67.Hc, 78.55.Cr, 68.35.Ct, 78.47.jd 



I. INTRODUCTION 



The kinetics of exciton recombination in semiconduc- 
tor QDs is often analyzed in terms of an exponential de- 
cay with one characteristic recombination time»i How- 
ever, the luminescence decay in QDs is typically non- 
exponentiall^'^, for which there are several reasons such 
as the contribution of dark excitons to the emission or 
the influence of Coulomb correlation effects. Neverthe- 
less, for a single QD in the strong confinement regime 
the bright exciton PL is found to decay monoexponen- 
tially.^''* For an ensemble of such QDs, on the other hand, 
nonexponential decays are often found, and a statistical 
analysis of the time-resolved emission demonstrates that 
this behavior can be attributed to a dispersion of radia- 
tive and/or nonradiative lifetimes of QD confined exci- 
tons. ■^'^ This ensemble decay at a specific energy results 
from the superposition of monoexponential PL decays of 
excitons which are localized in QDs having different sizes, 
shapes, and compositions.^ In case of continuously dis- 
tributed lifetimes t of excitons, which are characterized 
by the same recombination energy, their PL decay can 
be described by a distribution function G{t). 

We demonstrated recently that a nonexponential long- 
time decay of the exciton PL is characteristic of indirect 
band-gap (In, Al) As/ AlAs QDs with type-I band align- 
ment j^i^ In these structures the conduction band mini- 
mum is around the X- valley, while the valence band max- 
imum is around the F-point. In QDs the momentum is 
no longer a good quantum number, but the wave func- 
tion is distributed in momentum space over a range of 
k-vectors that is inversely proportional to the quantum 
dot size. This extension is, however, still much smaller 
than the separation between the F- and the X-valley, so 
that the indirect character of the band gap is maintained. 



As a result, direct band-to-band transitions of electrons 
resulting in an emission of a photon are strongly sup- 
pressed. Instead, the radiative recombination requires 
the involvement of scattering by phonons or at the het- 
erointerface, as has been demonstrated in indirect band- 
gap GaAs/AlAs and InAs/AlAs quantum wellSf^iiS For 
(In,Al)As/AlAs QDs it has been ascertained that the ra- 
diative exciton recombination is mainly caused by the 
scattering at the heterointerface between the (In,Al)As 
QD and AlAs matrixi^iii^ Hence, the exciton recombi- 
nation dynamics, namely the recombination time r and 
the lifetime distribution G{t) can yield valuable informa- 
tion on this interface. 

In this paper the dynamics of the exciton recombina- 
tion in ensembles of indirect band-gap (In,Al)As/AlAs 
QDs with varying sharpness of the QD/matrix inter- 
face is studied by time-resolved PL. We demonstrate 
that the radiative lifetime of the indirect in momentum- 
space exciton is strongly infiuenced by this sharpness. 
The decay can be well described by a power-law function 
I{t) ^ (1/i)", which can be accounted by a phenomeno- 
logical distribution function G{t) based on a single fitting 
parameter. 



II. SAMPLES AND EXPERIMENT 

The studied self-assembled (In,Al)As QDs, embed- 
ded in an AlAs matrix, were grown by molecular-beam 
epitaxy (Riber-32P system) on semi-insulating (001)- 
oriented GaAs substrates. The structures have one QD 
sheet sandwiched between 50-nm-thick AlAs layers grown 
on top of a 2G0-nm-thick GaAs buffer layer. The nomi- 
nal amount of deposited In As was about 2.5 monolayers. 
A 20-nm-thick GaAs cap layer protects the AlAs layer 



against oxidation. 

Recently, we demonstrated that diameter, density, and 
composition of (In,Al)As/AlAs QDs are determined by 
the growth conditions such as substrate temperature Tg, 
and growth interruption time ici'^ Three structures SI, 
S2, and S3 were grown for this study using the conditions 
hsted in Table HI According to these conditions the struc- 
tures have different (In,Al)As QD alloy compositions^^ 
However, as demonstrated repeatedly, despite of the in- 
termixing on the QD composition during the epitaxy, 
as-grown self-assembled QDs have a sharp QD/matrix 
interfacej ^"^!^^ The reason for the sharp interface forma- 
tion arises from the Stranski-Krastanov growth mode. 
The QD composition is determined by the intermixing 
during the dot formation due to mass transfer along the 
wetting layer J^"— The interfaces of QDs being indepen- 
dent on their composition are given by stable crystal- 
lographic planes, which ensure minimization of the in- 
terface energy. These planes provide interface stability 
against the intermixing with the matrix material during 
the overgrowth of the QDs. Nevertheless, the sharpness 
of the (In,Al)As/AlAs interfaces can be smeared out by 
means of high-temperature post-growth thermal anneal- 
ingpi^ Two of the studied structures, S2 and S3, were 
annealed for one minute at elevated temperatures Tan- 
Data on Tan are listed in Table U whereas technical de- 
tails can be found in Ref. [l^ The following, we will refer 
the annealed S2 structure as the S4 structure. 

The QD size and density are studied by transmission 
electron microscopy (TEM) using a JEM-4000EX sys- 
tem operated at 250 keV acceleration voltage. From the 
TEM-images we find that the self- assembled (In,Al)As 
QDs are lens-shaped with a typical aspect ratio (height 
to diameter) of 1:44 TEM plane- view images and the re- 
spective histograms of the QD diameter distribution are 
shown in Fig. [1] for all structures. The average diame- 
ters T'av and the diameters corresponding to the larger 
Dl and smaller Dg half-widths of the QD size distribu- 
tion obtained from the TEM images, are summarized in 
Table |T1 Additionally, the size dispersion Su is listed, 
which is defined by the ratio of the half-width of the 
Gaussian distribution of QD sizes to the average diame- 
ter 5*0 = 100% X {Dl - Ds)/Dav- 

The unannealed structure SI has a relatively narrow 
distribution (5'd = 40%) around a small QD size of 
T'av = 5.5 nm. The structures S2, S3, and S4 contain 
QDs with larger diameters, which are also distributed 
over broader ranges. The largest lateral QD sizes are 
found for structure S4, which was annealed at the high- 
est temperature of 800° C. The comparison between the 
annealed structure S4 and the unannealed structure S2 
indicates that post-growth annealing leads to an increase 
of the QD diameter with a broader distribution, in good 
agreement with previous results J4 The annealing results 
in an increase of the average QD diameter from 13.8 to 
19.6 nm, and 5d from 60% to 75%, respectively. The 
diameter increase by annealing is a result of InAs dif- 
fusion from the QD into the surrounding AlAs matrix. 
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FIG. 1. (Color online) TEM plane-view images (upper panels) 
and histograms of the QD-diameter distribution (lower pan- 
els), including the size dispersion fitted by Gaussian curves 
for the structures S1-S4. The histograms show the incidences 
of finding a QD with a specific diameter in a TEM image 
within an area ranging from 0.3 to 0.5 fira^ . The half- widths 
of the QD-size distribution are marked by the horizontal lines. 
The average diameter Dav and the large Dl and small Ds 
half-widths of the QD-size distribution are indicated for the 
structure S3. For the structure S2, _D£ is taken as the diam- 
eter corresponding to indirect band gap QDs, see caption to 
the Table III 



Therefore, the annealing results in the appearance of a 
diffused Ina;Ali_a;As layer around the QD/matrix inter- 
face. It is obvious, that the thickness of this In^Ali_^As 
layer depends on the annealing temperature and dura- 
tion. The sharpness of the interface is defined as the de- 
gree of the spatial separation between the different ma- 
terials of (In,Al)As and AlAs. Thus, the sharpness of 
the QD/matrix interface can be described in terms of 
the thickness of this diffused layer: a sharp (blurred) in- 



TABLE I. Growth parameters and annealing temperatures for the studied (In,Al)As/AlAs QDs. The average diameter (Dav), 
diameters corresponding to the smaller (-Ds) and larger (-Dl) half-width of the QD size distribution, the QD density, the size 
dispersion (So), and the composition of the (In,Al)As/AlAs QDs are given as well. Additionally, the exponent of the PL 
decay curve a, and the parameter to in the exciton lifetime distribution described by Eq.([3| are listed. Note, the relation 
7 = + 1 defines another parameter of the exciton lifetime distribution 7. For the structure S2 the diameter corresponding to 
the larger half-maximum in the QD size distribution is related to QDs with direct band gap. In order to have for structure S2 
an additional ensemble of indirect band gap QDs with a characteristic diameter being different from Dav and Ds , we choose 
the diameter (marked in Fig. [T]by D£ = 17 nm) that is related to indirect band gap QDs. The parameter values marked in 
the Table with an asterisk belong to the QD ensemble with the characteristic diameter D£. 
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terface corresponds to a thin (thick) In^^Ali^j^As layer. 
Therefore, the structures S1-S4 provide us with a repre- 
sentative set of QD ensembles having different diameters 
and interface sharpness. 

The steady-state and time-resolved PL measurements 
were performed at a temperature of T = 5 K. For the 
excitation of the steady-state PL a He-Cd laser with a 
photon energy of 3.81 eV was used. The time-resolved 
PL experiments were established by the third harmonic 
of a Q-switched Nd:YV04 laser (3.49 eV) with a pulse du- 
ration of 5 ns. The pulse-repetition frequency was varied 
from 1 to 100 kHz and the pulse energy density was cho- 
sen between 40 nJ/cm^ to 12 /iJ/cm^. The emitted light 
was dispersed by a 0.5 m monochromator and detected by 
a GaAs photomultiplier operating in the time-correlated 
photon-counting mode. In order to monitor the PL decay 
in a wide temporal range up to 0.5 ms the time resolu- 
tion of the detection system was varied between 1.6 ns 
and 200 ns. 



III. EXPERIMENTAL RESULTS 

A. Steady-state photoluminescence 

Normalized PL spectra of the studied structures are 
shown in Fig. [2] Two PL bands, denoted by DQD 
and IQD, are observed in the spectra. Recently, we 
have shown that the PL kinetics of the DQD and IQD 
bands are considerably different. The DQD and IQD 
bands result from exciton recombination in QDs with di- 
rect and indirect band gap, respectively. Hereby, it has 
been evidenced by pulsed excitation measurements that 
the intensity of the DQD emission band drops quickly 
with a decay time shorter than 20 ns, while the IQD 
band decays nonexponentially during hundreds of /Lts4 
These two types of QDs coexist in these ensembles of 
(In,Al)As/AlAs QDsi'-ii For the PL spectra of the SI, 



S2, and S3 structures the direct DQD band has been 
observed at energies below 1.65 eV. This energy well co- 
incides with the observed and calculated boundary be- 
tween the direct and indirect band-gap (In,Al)As/AlAs 
QDs.^^ 

Since the shape of the PL emission reflects the distri- 
bution of QD sizes, we establish in the following the re- 
lation between the parameters characterizing the spectra 
and the geometrical quantities of the average diameter 
-Dav and size dispersion Sd. This relation is given by 
the following features: 

• The increase of Dav and 5d for the as-grown struc- 
tures SI and S2 leads to a low-energy shift from 
1.8 to 1.7 eV and an increase in full width at half 
maximum (FWHM) of the IQD band from 115 to 
195 meV, respectively. Additionally, the intensity 
of the DQD band is increased by one order of mag- 
nitude in S2 compared to SI. 

• Increase in Dj^y and Sd, as a result of high temper- 
ature annealing of the S2 structure and its transfor- 
mation to the S4 structure, leads for IQD band to 
high-energy shift from 1.7 to 1.85 eV and decrease 
of FWHM from 195 to 75 meV. The annealing re- 
sults also in decrease of the DQD band intensity in 
the S3 structure and disappearance of such band in 
the S4 structure. 

In order to explain the obtained results we have to take 
into account that the energy of the optical transition is 
determined by two factors: (i) the quantum confinement 
energy which decreases with increasing QD size; and (ii) 
the band-gap energy of the (In,Al)As alloy in the QD 
which increases with decreasing InAs fraction. In lens- 
shaped QDs the confinement energy is mainly determined 
by the QD height. The average QD composition can be 
evaluated from comparison of the IQD band energy po- 
sition with results from model calculations.^^ The deter- 
mined QD compositions are collected in Table HI 
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FIG. 2. (Color online) Normalized PL spectra of the different 
(In,Al)As/AlAs QD structures SI to S4, excited by a He-Cd 
laser with a power density of 10 W/cm^. The PL intensity 
of the SI structure is multiplied by a factor 0.7 for better 
visibility. The vertical dotted area marks the boundary be- 
tween the emission from QD excitons either direct or indirect 
in momentum-space as calculated in Refs.[Tll, @. The vertical 
arrow indicates the energy of exciton recombination in QDs 
with typical diameter _D£ for the structure S2 (see caption 
of Table |I|. Schemes on top: Real space band structures of 
direct and indirect band gap QDs, including the energetically 
lowest F- and X-conduction band states as well as the heavy- 
hole (hh) states. 



The comparison of the observed optical transition en- 
ergy with calculated ones shows that the low-energy shift 
and broadening of the IQD band, going from the SI to the 
S2 structure, are caused by a decrease in quantum con- 
finement and an increase in the size dispersion 5*0, while 
the change in the QD composition is negligible. On the 
other hand, the high-energy shift of the IQD band for 
the S4 structure compared to S2 is due to an increase 
of the (In,Al)As alloy band-gap energy with decreasing 
InAs fraction in the QD alloy composition. This compen- 
sates the decrease in quantum confinement energy due to 
the annealing induced increase of the QDs height based 
on a fixed aspect ratio. The at first sight unusual reduc- 
tion of the FWHM of the IQD band with increasing ^d 
was explained in our previous studyJ^ 

We calculated in Refs. [ill, |^ the energy separation be- 
tween the optical transitions of direct and indirect exci- 
tons and showed that it weakly depends on the QD size, 
shape, and composition. Therefore, the change in the 
relative intensity of the DQD band reflects the change 
in the relative layer density of the direct band gap QDs. 
This density is a function of the QD size and composition. 
The disappearance of the DQD band in the emission of 
the S4 structure results from the lower InAs fraction in 
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FIG. 3. (Color online) Low temperature (T = 5 K) PL ki- 
netics obtained at the maximum of the IQD band (1.856 eV) 
for the S4 structure using different energy densities of excita- 
tion pulse (note the double- logarithmic scale). The excitation 
pulse ends at t = 10 ns. The pulse repetition frequency is 
1.5 kHz, which is sufficiently low to monitor a PL intensity 
decrease by five orders of magnitude between sequent laser 
pulses. From top to bottom the excitation pulse densities P 
given in nj/cm^ are: 1.2 x lO", 4 x 10^, 1.2 x 10^, 400, 120, 
and 40. The inset demonstrates the PL kinetics measured 
at P = 900 nJ/cm^ for different pulse repetition frequencies, 
from top to bottom in kHz: 100, 30, 10, and 1. While for the 
power density a logarithmic scale is again chosen, the inten- 
sity axis is scaled linearly, in order to underline the changes 
in the PL decay. 



the QD alloy composition. This, in turn, gives rise to a 
conversion of the band gap from a direct to an indirect 
one. 



B. Time-resolved photoluminescence 

For direct band gap (In,Ga)As QDs with typical ex- 
citon lifetimes of about 1 ns the condition that the QD 
exciton population does not exceed one exciton can be 
easily established. As an example, for (In,Ga)As/GaAs 
QD ensembles, which are excited by picosecond pulses (at 
13.2 ns pulse separation), the average exciton population 
per dot is smaller than 0.15, when an average excita- 
tion density of 8 W/cm^ (energy density 100 nJ/cm^ per 
pulse) is used.^ 

However, in indirect band gap QDs with long exciton 
lifetimes the optical excitation has to be carefully cho- 
sen in order to avoid accumulation of electron-hole pairs 
and formation of multiexciton complexes in the QDs. 
For that purpose, specific experimental conditions need 
to be established. The number of excitons, which are 
photogenerated in the QD surrounding matrix and cap- 
tured in the QDs per laser pulse, is mainly determined 
by the pulse energy density, but is independent of the 
QD band gap structure because the relaxation from the 
excited F to the X ground state is very fast. The repeti- 






(b) ' 


i..,.^^^ low 


S2 - 


_rrax N^ 




li]|^|[^\^\^^ 




"''°'hi^^\\/^ 


V\ ■ 


1 1 1 


\v\ 



-! 




low 


(d) 
S4 - 


xi6 


^ 


" x1 


max 


xO.2 


high 








1 



10 10 10 
Time (ns) 



10' 



10 10 10 
Time (ns) 



10= 



FIG. 4. (Color online) Low temperature (T = 5 K) PL kinet- 
ics measured for the structures S1-S4 at the intensity maxi- 
mum ('max') and at half of this maximum on either the high 
energy ('high') or the low energy ('low') side of the IQD emis- 
sion band. The excitation density amounts to P = 40 nj/cm^ 
and the pulse repetition frequency is 1.5 kHz. The excitation 
pulse ends at 10 ns. The scaling factors of the PL intensities 
are introduced for better visualization. The thick solid lines 
show the modeled results with the distribution of G{t) de- 
scribed by Eq. (3) with the parameters presented in Table |T| 
For the S2 structure curve 'low' corresponds to the exciton 
recombination in QDs with typical diameter _D£. 



tion frequency of the excitation pulses should be reduced 
to a level that there is sufficient time for the excitons 
to recombine between subsequent pulses, so that multi- 
exciton complexes are not created. Since the lifetime of 
the indirect excitons in the (In,Al)As/AlAs QDs exceeds 
the one of direct excitons in (In,Ga)As/GaAs QDs by 
up to five orders of magnitude?., the pulse repetition fre- 
quency as well as the average excitation density should 
be decreased correspondingly. To ensure that we actually 
study the recombination dynamics of single excitons in 
the (In,Al)As/AlAs QDs, the PL kinetics was measured 
at different excitation densities and repetition frequen- 
cies. 

Figure |3| shows low temperature PL kinetics measured 
at the IQD band maximum of the S4 structure for dif- 
ferent excitation pulse energy densities. The kinetics of 
the other structures is similar to the presented one. The 
transient PL data are plotted on a double logarithmic 
scale, which is convenient to illustrate the nonexponen- 
tial character of the decay over a wide range of times and 
PL intensities. The recombination kinetics demonstrates 
two distinctive stages: (i) a relatively flat PL decay im- 
mediately after the excitation pulse up to approximately 
1 fis and, subsequently, (ii) a reduction in the PL in- 



tensity which can be described by a power-law function 
I{t) ^ (1/0" ^s shown in our previous studies^ One can 
see that in the case of the stage (i) , a high power excita- 
tion results in a fast decay of the exciton PL. It can be as- 
signed to the recombination of multi-exciton complexes)^ 
By decreasing the power down to P = 120 nJ/cm-^ the 
decay decelerates, below this power the intensity does 
not temporally change thus indicating a saturation level. 
It is induced by the recombination of single excitons in 
the QDs. Taking into account the absorption of the laser 
light in the AlAs matri x^ -' i-" ' together with the QD den- 
sity, the average number of excitons captured in a QD 
per pulse is estimated to about 0.3 for P = 120 nJ/cm^. 
An increase in the repetition rate of the excitation pulses 
at fixed pulse power also results in an acceleration of 
the initial kinetics stage, as shown in the inset of Fig. [31 
As the sample has to be excited by each pulse when it 
has reached its equilibrium state, for the further studies 
presented in this paper we select P = 40 nJ/cm^ . It cor- 
responds to an average QD exciton-population of 0.1 per 
pulse at a repetition frequency of 1.5 kHz which equals 
to a time interval of 670 fis between the sequent pulses. 

Fig. m shows PL decays for the structures S1-S4 which 
have been measured at different energies. The selection 
of different detection energies provides information on 
the exciton recombination in QDs with different charac- 
teristic sizes in the ensemble. Hereby, the energy of the 
intensity maximum of the IQD band (the curves labelled 
'max' corresponding to recombination in QDs with the 
diameter -Dav) and of the half maximum (curves 'high' 
and 'low' related to recombination in QDs with diame- 
ters Ds and I?l, respectively). The following features 
of the PL decay can be extracted: (i) for each structure 
the exponent a of the power-law decay is determined by 
fitting of the second stage of the decay curves with the 
form I{t) r^ (1/t)", as listed in Table ID The exponent 
a increases monotonically across the IQD band from the 
low to the high energy side, (ii) The second decay stage 
starts hundred nanoseconds after the end of the excita- 
tion pulse in the unannealed structures SI and S2, while 
for the annealed structures S3 and S4 it begins several 
microseconds after the pulse. The PL decay obviously 
depends on the sample characteristics and as it will be 
shown in the following, it is affected by both the typical 
QD size in the ensemble and the interface sharpness. 

To obtain a quantitative description of the effect of QD 
size and interface sharpness on the exciton lifetime in the 
(In,Al)As/AlAs QD ensembles, we need to construct the 
distribution function G{t), which controls the observed 
PL kinetics. 



C. Exciton-lifetime distribution in (In,Al)As/AlAs 
QD ensembles 

Nonexponential decays of the exciton PL intensity I{t) 
are frequently described by stretched-exponentials of the 
form I{t) ex exp[— (i/r)'^], including a constant lifetime 



T and a dispersion factor /3iS^i2^ The stretch parameter 
< /3 < 1 quaUtatively expresses the underlying dis- 
tribution function G{t): a broad distribution results in 
(3 <^ 1, while for a narrow one /3 is about 1. However, 
the evaluation of the lifetime distribution on the basis of 
the stretched-exponential model is mathematically com- 
plicated and feasible only for specific /3 values (see Ref. y 
and references therein). Alternatively, the distribution 
G{t) can be determined using the following equation: 



lit) 



I GiT)exp(--jdt. 



(1) 



Here, G{t) is established via either the numerical so- 
lution of the integral equation ([TJ,^^"^^ or an assumed 
analytical expression of G{t) with a set of fitting param- 
eters. Van Driel et ali^ assumed recently that the most 
successful distribution function to model G{t) in a QD 
ensemble among different analytical expressions like nor- 
mal, Lorentzian, etc distributions is a log-normal func- 
tion given by: 



G{t) = — exp -■ 



w 



(2) 



with the constant A and the maximum tq of the exci- 
ton lifetime distribution. The dimensionless parameter 
w describes the distribution width Aj^ /^ of the inverse 
recombination times at the 1/e level: Ai/^ — ^ sinh(it;). 
This distribution was successfully used to describe the 
nonexponential decay of the exciton PL intensity over 2-3 
orders of magnitude for different QD systems with a con- 
tinuous distribution of direct exciton lifetimes. Among 
them are ensembles of CdSe/ZnSe colloidal QDs^^ and 
dye molecules embedded in a photonic crystal. ^^ 

Since in (In,Al)As/AlAs QDs the excitons recombine 
via radiative recombination only^^ the nonexponential 
decay is the result of the dispersion of the excitonic ra- 
diative times in the ensemble. In order to determine G(t) 
we follow the approach of Ref. ^ and fit the PL kinetics 
with Eq. ([1]) using the log- normal distribution of Eq. © . 
Unfortunately, the log-normal distribution does not al- 
low us to describe the PL decay satisfactorily over the 
whole dynamical range. One can see in Fig. [5] for the 
structures SI and S4 that for the decay curves, measured 
at maximum intensity of the IQD band (corresponding 
to recombination in QDs with diameter I?av), this dis- 
tribution allows to fit either the initial stage (curves '1') 
or the long-time stage (curves '2') of the decay, using dif- 
ferent sets of parameters which are given in the figure 
caption. 

In order to describe the exciton PL decay in our struc- 
tures over the whole dynamical range of five orders of 
magnitude we propose a non-symmetric phenomenologi- 
cal distribution G{t), which is suitable for fitting power- 
law decays I{t) - {l/t)°': 
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FIG. 5. (Color online) PL decay of the samples SI and S4 
of the (In,Al) As/AlAs QDs measured at the IQD band maxi- 
mum for T = 5 K, the experimental data are shown by circles. 
The excitation density is P = 40 nj/cm^, and the pulse rep- 
etition frequency is 1.5 kHz. The end of the excitation pulse 
corresponds to the time 10 ns. The modeling of the PL de- 
cays by different distribution functions G(t) is shown by the 
lines. The corresponding distribution functions are given in 
the insets, (a) Structure SI: curves T' and '2' are based on 
model calculations using the log-normal distribution G{t) of 
Eq. ^ with parameter sets (ro = 0.2 fis, w = 1.55) and 
{to = 0.05 /is, w = 3.50), respectively. Curve '3' is for G{t) 
of Eq. (|3J with the parameters ro =0.1 /is, and 7 = 2.75. (b) 
Structure S4: curves T' and '2' are also modeled by G{t) of 
Eq. ((2} with the parameter sets (ro — 5.8 us, w = 1.35) and 
(ro = 3.0 /is, w = 2.10), respectively. Curve '3' is given by 
G{t) with ro = 5.2 /is, and 7 — 3.40. 



Here C is a constant, and tq characterizes the maximum 
of the distribution of exciton lifetimes. The parameter 
7 in the above Eq. ^ is defined as a -t- 1. By use of 
a double-logarithmic scale the power-law decay (1/i)" 
represents a line with slope a. Note, in knowledge of the 
exponent a of the experimental decay curve, only a single 
free parameter tq is required to describe the PL decay by 
Eq. ®. 

We use Eq. ^ to fit the decay curves in Fig. \5[ The 



experimentally determined a values (see Table H]) yield 
the values of 7 — 2.75 and 3.40 for the structures SI 
and S4, respectively. One can see excellent agreement 
between the experimental data and the calculations for 
the whole PL decay comprising five orders of magnitude 
from the low nanosecond to the high microsecond region, 
shown by the curves '3' in Figs. [SJa) and [SJb). The 
fits were obtained for the G{t) of Eq. (|3]), depicted in 
the insets of Fig. [5] by the curves labeled '3'. The used 
fit parameters are tq — 0.1 /is and tq = 5.2 /is for the 
structures SI and S4. 

Let us discuss now the parameter 7 in Eq. ([3]). The 
comparison of the distribution functions Eqs. ([2]) and ([3|), 
as shown by the curves '1' and '3' in the insets to Fig. [SJ 
clearly demonstrates that the strong difference in G(t) 
for lifetimes smaller than tq has very little influence on 
the initial stage of the PL decay (tq > t). Thus, the 
decay curves are mainly contributed by recombination of 
excitons with lifetimes exceeding tq. Therefore, the value 
of the parameter 7, which is the exponent of the long- 
lifetime tail of G(t)^ can be used as qualitative measure 
for the effective width of the distribution Gir). Accord- 
ing to Eq. ([3]), an increase in 7 reduces the dispersion of 
the r values which contribute to the long-time tail of the 
kinetics curve. 

We would also like to demonstrate that the filling of 
the QDs with multi-exciton complexes at high excita- 
tion densities distorts G(r). The Fig. [5] illustrates the 
G(t) distributions obtained via fitting of the recombi- 
nation kinetics for the exciton population of less than 
one (40 nJ/cm^), and in case of multi-exciton occupation 
(1 /iJ/cm^). One can clearly see in Fig. IH] that multi- 
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FIG. 6. (Color online) Normalized distribution functions G(t) 
corresponding to QD ensembles with the same average diam- 
eter 73 Av = 19.8 nm (structure S4). The two dependences are 
obtained by fitting the recombination kinetics measured for 
different excitation energy densities: 40 nj/cm^ (black solid 
line) and 1 /iJ/cm^ (red dashed line). 



exciton QD occupation shifts the distribution maximum 
to shorter times. Thus, in order to reveal the G{t'), be- 
ing intrinsic for the ensemble of indirect band gap QDs, 
the average exciton population in the dots should be less 
than unity. 



D. Effect of interface sharpness on exciton lifetime 

Using Eqs. ([T]) and Q we determine the G(t) distribu- 
tions for the four studied structures by fitting the decay 
curves presented in Fig. 21 The values of the parameters 
7 = a -|- 1 , and ro resulting from the fitting are collected 
in Table HI The results of the fitting are shown in Fig. |4] 
by solid lines. 

The following conclusions on the distribution proper- 
ties can be drawn from the data in Table ID 

• A monotonic decrease of tq and 7 with increasing 
QD diameter (from Dg to D\) is a common feature 
for each structure. 

• The relative change of tq and 7 with changing QD 
diameter from Da, to Z3l is larger in the as-grown 
structures SI and S2 than in the annealed struc- 
tures S3 and S4. This is also evidenced by Fig. [7] (a) 
for the structures S2 (with smaller size dispersion 
Sd — 60%) and S4 (with larger size dispersion 
Su = 75%). 

• Despite of the larger QD diameters in the annealed 
structures S3 and S4 than in the as-grown struc- 
tures SI and S2, tq is much smaller in SI and S2 
than in S3 and S4. Also, it is independent of the 
QD diameter, as demonstrated in Fig. [7] (b) for the 
QDs with their average diameters 13 av in structure 
SI to S4. 

These features give us the possibility to distinguish be- 
tween the effect of QD size and interface sharpness on 
the exciton lifetime. One can see in Fig. [7] (b) and in Ta- 
ble |T] that, despite of the large difference in QD diameter 
(Dav — 5.5 and 13.8 nm for SI and S2), these struc- 
tures have similar distributions of exciton lifetime. On 
the other hand, the comparison of G(t) in QD ensembles 
with similar diameters (Dl= 17 nm, Z^av = 18.3 nm, 
and Dav = 19-6 nm for the structures S2, S3, and S4, 
respectively) highlights a strong increase of the exciton 
lifetime with decreasing interface sharpness by about two 
orders of magnitude. Thus, the exciton lifetime in indi- 
rect band gap (In,Al)As/AlAs QDs is mainly determined 
by the interface sharpness, while its dependence on the 
QD size is much weaker. Nevertheless, the monotonic 
decrease of tq with increasing QD diameter from Dg to 
-Dl for each of the studied structures indicates that effect 
of QDs size on recombination time is also important. 

Our theoretical estimations are in good accordance 
with the experimental data. As demonstrated in the Ap- 
pendix, the lifetime of an indirect in momentum space 
exciton can be described by the form t oc expid / a -^ d / L) 



8 



due to momentum scattering at the interface. Here, 

a is the lattice constant, L is the QD height, and d 

is the thickness of the diffused (In,Al)As layer at the 
QD/matrix interface. 

Note, the ratio d/a > 1 increases with the increase of 
the thickness d of the diffused layer. Therefore, the exci- 
ton lifetime is mainly determined by d/a. It is reasonable 
to assume that d/a changes weakly with the QD size for 
a particular structure type. Hence, the dependence of 
the exciton lifetime for such a structure type on the QD 
size is specified by the second term d/L < 1. As a re- 
sult, the exciton lifetime decreases in agreement with the 
experimental data, for example for the S2 structure in 
Fig. [3 (a). 
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FIG. 7. (Color online) Distributions G{t) corresponding to 
(a) the QDs with diameters Dl = 17 nm, Dav = 13.8 nm, 
and Ds = 9 nm for structure S2 (from left to right); and the 
QDs with typical diameters Dl = 28 nm, -Dav = 19.8 nm, 
and Ds = 12 nm for structure S4; (b) the QDs with Dav 
diameter 5.5 nm in structure SI, 13.8 nm in structure S2, 
18.3 nm in structure S3, and 19.6 nm in structure S4. 



Decrease of exciton lifetime with increase in d/a value 
is restricted by a rate of phonon emission. When the 
rate of non-phonon recombination becomes smaller than 
that of the phonon emission the exciton recombination 
is mainly determined by the phonon emission. Actu- 
ally, we demonstrate in Ref. |l8| that very high anneal- 
ing temperature (950° C), which leads to very smooth 
QD/matrix interface, results in appearance of phonon 
replicas in the PL spectrum of the annealed structures 
with In(Al)As/AlAs QDs. 

Let us now consider qualitatively the effect of interface 
sharpness on the effective width of the G{t) distribution. 
The exciton lifetime is proportional to exp(ci/a -I- d/L). 
Therefore, the width of the lifetime distribution should be 
determined by variations of the argument d/a+d/ L value 
for different QDs. The dispersion of a varying quantity 
is inversely proportional to the square root of its mean 
valuCi^ Thus, the width of G{t) should decrease with in- 
creasing value of d/a -\- d/L. The experimental data tend 
to confirm this expectation. Indeed, the effective width 
of the G{t) distribution (which is inversely proportional 
to 7) for QDs with similar sizes decreases with increas- 
ing thickness of the diffused layer. While an increase of 
the ratio d/L with increasing QD size for a fixed ratio 
d/a results in a reduction of d/a + d/L, and, thus, in an 
enhancement of the G{t) width (decrease of the value 7 
is shown in Table H]) . 



IV. CONCLUSION 

The dynamics of the exciton radiative recombination in 
(In,Al)As/AlAs QD ensembles with a type-I band align- 
ment, but an indirect band gap has been investigated. 
Due to the different dot sizes and/or QD/matrix inter- 
face sharpness the exciton recombination dynamics shows 
a nonexponential decay behavior that can be described 
by a power- law function as a result of the superposition of 
multiple monoexponential PL decays with different life- 
times. The lifetime of these excitons is mostly deter- 
mined by the thickness of the diffused (In,Al)As layer 
at the QD/matrix interface, while its dependence on the 
QD size is weaker. We have proposed a phenomenological 
equation for the distribution G{t) of the radiative exci- 
ton recombination times in such QD ensembles, which 
can describe the power-law PL decay over five orders of 
magnitude very well, using one fitting parameter only. 
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APPENDIX 



The exciton wave function can be written as 



(re,rh) = /e(re)//i(rh)J(re - rh) 



(4) 



where rg and rh are the coordinates of electron and 
hole, fe and //,. are their wave functions in absence of 
the electron-hole Coulomb interaction, while J(re — rh) 
takes into account this interaction. The exciton recom- 
bination rate is proportional to ^^(0)|(/eV//i)p. If we 
decompose the wave functions into products of the Bloch 
waves times envelope wave functions Fe{r) and F/i(r) and 
assume that V acts only on the Bloch amplitudes, then 
we find that the exciton recombination rate is propor- 
tional to the square of module of the convolution of the 
envelopes: 



(5) 



r= / Fe(r)F,:(r)d^r. 



To estimate this integral we assume the following form 
for the envelopes in the vicinity of the interface between 
QD and matrix (z = 0): 

Fe{z) = exp{iqz) yj A{k) exp(ifcz), (6) 

fc 

F„(z)-;^B(p)exp(zpz), (7) 

p 

where q — T:/a, and A, B are coefficients, which can be 
obtained from the boundary conditions. They are con- 
stants in the infinite crystal, but depend on the electron 
and hole momenta k and p in the quantum dots. Sum- 
ming on these values spreads over k, p ^ 1/L where 
L ^ a is the dot size (for our lens-shaped QDs L is the 
height of QD). Thus, we can assume k, p <^ q and the 
value of the exciton recombination rate at the QD/matrix 
interface is determined by the integral: 



Integration over the dot interface leads to estimation of 
r^ as the ratio of the number of atoms located at the 
interface to the total number of atoms in the quantum 
dot. Estimation of F^ at a diffused interface can be done 
assuming that p and k vary smoothly with position inside 
the interface layer: p(z) = ^2rah{E — U{z)) and k{z) — 
where U{z) is the potential profile of 
Then 



^27n,{E~U{z)), 
a diffused interface. 



r. = 



+ CXD 



exp(z(j^(z) — k{z) + q)z)dz. 



(10) 



To estimate F^ we consider the integrand in Eq. (jlOl) as 
a function of the complex variable z. We can displace 
the integration contour from the real axis z into the up- 
per half-plane, up to the nearest singularity Zp of the 
potential U{z). This value is about id, where d is the 
characteristic thickness of the diffused interface. For 



U{z) 



Uo 



1 + cxp(— z/d) 



(11) 



the value of Zp is given by iTrd; the actual value of Zp 
depends on the model for the interface. Therefore, the 
integral (I10[) contains the exponential factor exp{iqzp) ~ 
exTp{qd). Evaluation of the integral similar as done for 
Eq. Uni) in Ref. i32 results in 



exp(— gd) a 

^ — exp 

q TT \ a 



(12) 



Note that estimations ([8))- p0)) suppose the large size of 
the QD. Decrease of the QD size, e.g. its height, leads to 
increase of the electron and decrease of the hole energies 
due to the size quantization as well as non-zero value of 
the electron and hole momenta (which is about ■nh/L). 
To take this fact into account, we have to add ii/L to the 
exponent Eq. (|10l) . i.e. replace q with q + ir/L in Eq. (flT 
and Eq. ^^. Then Eq. dH]) accepts the form 



exp[i(p — fc + q)z\dz, 



(8) 



which is zero far from the interface, because of the os- 
cillating factor exp(i(7z). We can estimate F^ at a sharp 
interface as 

T z = I exp[i(p — fc -I- q)z\dz ^ — ^ — . (9) 

J iq TT 



Tz ^ — exp 

TT 



■ii^-i) 



(13) 



Taken into account that exciton recombination rate is 
inversely proportional to exciton lifetime our estimation 
demonstrates that r '-^ exp [^ -I- -^j i.e. increase in thick- 
ness of diffused layer at the QD/matrix interface results 
in exponential increase of exciton lifetime. 
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